The Dynamical Dipole Mode in Dissipative Heavy Ion Collisions 
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I We study the effect of a direct Giant Dipole Resonance (GDR) excitation 

Q : 

in intermediate dinuclear systems with exotic shape and charge distributions 
formed in charge asymmetric fusion entrance channels. A related enhance- 
ment of the GDR gamma yield in the evaporation cascade of the fused nucleus 
is expected. The dynamical origin of such GDR extra strength will show up in 

X' 

S-h ' a characteristic anisotropy of the dipole gamma-emission. A fully microscopic 

analysis of the fusion dynamics is performed with quantitative predictions of 
the GDR photon yield based on a dynamics- statistics coupling model. In 
particular we focus our attention on the energy and mass dependence of the 
effect. We suggest a series of new experiments, in particular some optimal 
entrance channel conditions. We stress the importance of using the new avail- 
able radioactive beams. 
PACS: 24.30.Cz, 25.70.Jj, 25.70.Lm, 25.60.Pj 
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I. INTRODUCTION 



The isovector giant dipole resonance represents a well established collective motion of 
finite nuclei extensively studied since more than fifty years ffl. The dependence on the nu- 
clear structure of the reference state on top of which the collective mode is built has already 
suggested the use of GDR properties in order to study nuclei far from the ground state. 
The time structure of the GDR mode actually allows a possibility of using it as a probe of 
nuclear systems very far from normal conditions, see the nice reviews 0J3|. 
We can roughly estimate an oscillation period of ^h- ~ 80 - 100/m/c and a mean life-time 
Ti/Tgdr — 50/m/c. Since the spreading width is the largest contribution to the total width 
^gdr, 50/m/c is roughly also the time needed to build up the collective GDR mode in a 
Compound Nucleus (C.N.). All these time scales are relatively short and this makes the 
GDR an ideal probe to study nuclear systems under extreme conditions. 
Charge equilibration takes place on time scales of few units in 10 2 fm/c: therefore in charge 
asymmetric entrance channels we expect to see a "direct dipole" collective excitation, with 
special dynamical features, which should give an extra yield to the photon spectrum in the 
GDR region. The idea is that we can form a fused dinuclear system with the charge not yet 
equilibrated and therefore with some extra dipole strength of non statistical nature. In this 
sense the " Dynamical Dipole" is a direct dipole contribution present in dissipative collisions 
initiated in charge asymmetric entrance channels. This new pre-equilibrium collective mode 
has been predicted some years ago @,|] and recently clearly observed in fusion |||| and 
deep-inelastic PJlOll reactions. Some experiments have been performed also at intermedi- 



ate energies, where other dynamical entrance channel features can be present [[Lljjl^l . The 
appearance of this prompt dipole effect gives important information on the charge equili- 
bration dynamics in connection to the reaction mechanism. A fully microscopic analysis 
is performed with quantitative predictions of the GDR photon yield based on a dynamics- 
statistics coupling model. Independent information on the damping of the dipole mode in 
very excited nuclear systems can be derived. Moreover the coupling between isovector dipole 
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and large amplitude isoscalar monopole oscillations can be studied, of large interest also for 
the discussion on the fate of the GDR in hot nuclei. Since the effect is strongly dependent 
on fusion dynamics we study the beam energy variation of this extra dipole strength and we 
suggest some optimal conditions for the relative observation. 

II. GDR IN CHARGE ASYMMETRIC FUSION: THE DYNAMICAL DIPOLE 

Fusion processes in charge asymmetric entrance channels give a clear example on how 
emitted GDR photons carry information on a Compound Nucleus (CN), formed during the 
reaction dynamics, very far from normal conditions. 

In dissipative collisions energy and angular momentum are quickly distributed among all 
single particle degrees of freedom while charge equilibration takes place on larger time scales 
fl"3|j . Therefore for charge asymmetric entrance channels at the time of CN formation we 
can easily have some relic of a pre-equilibrium GDR from dipole oscillations in the isospin 
transfer dynamics |J. The result will be an enhanced GDR photon emission with special 
features due to the non-statistical nature of this extra contribution. 

An enhancement of the GDR gamma ray yield has been observed experimentally in 
fusion reaction between heavy ions with different neutron to proton ratios. Flibotte et al. 
@, in collision between 40 C a (N/Z = 1) and 100 Mo, (N/Z = 1.38), at AAMeV beam energy, 
producing a CN at 70MeV excitation energy, counted a number of emitted GDR photons 
over the whole cascade 16% larger than in the case when the same CN was formed with a 
symmetric N/Z combination. Similar effect has been observed by Cinausero et al. 0: an 
increase of 36% resulted when they compared GDR photon yields from a CN populated at 
110 MeV excitation energy via the reactions 16 O {N/Z = 1) + 98 Mo (N/Z = 1.33) and 50 Ti 
(N/Z = 1.18) + 64 Ni (N/Z = 1.28) respectively. 

In fact, theoretically, an extra-yield of gamma rays was predicted some years ago by 
Chomaz et al., ||, in a simple model of a GDR phonon gas interacting with the CN. The 
essential quantity in this approach is the mean number of excited GDR phonons at the time 
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when a CN is formed, riQ DR . The emission of GDR photons is enhanced if tvq DR is greater 
than the value corresponding to a statistical equilibrium between the GDR oscillator and 
the CN heat bath. The N/Z asymmetry of colliding heavy ions will certainly trigger in 
the early stage of the reaction the charge equilibration process. Several experiments and 
realistic simulations have indicated that the related dipole oscillation has the characteristics 
of a GDR type quantal collective motion \TM , WB , ITfl] , P7| , Expressed in terms of 



the Brink- Axel hypothesis |19j this means that a GDR may be built not only on equilibrium 
states, of a warm nucleus for example, but also on non-equilibrium states during the CN 
formation phase. n < Q > DR represents just the number of phonons in this mode at the time of 
the CN formation, tcN- 

Although it was assumed that the quantity tiq dr is intimately related to the presence of 
pre-equilibrium effects, the way it is determined and affected by the projectile-target N/Z 
asymmetry and entrance-channel dynamics has not been investigated yet. From the previous 
discussion we clearly distinguish three main phases in the entrance channel dynamics: I - 
The approaching phase, with the two partners still keeping their own response properties; 
77 - The dinuclear phase, with relative collective response; 777- The CN formation. We 
will call t = the starting time of the phase II, i.e. the onset of the main dissipation 
mechanisms, including the charge equilibration. In this phase a pre-equilibrium dipole mode 
is present in the composite system, damped with a spreading width ]> = hfj,(t), where the 
time dependence is expressing a non-equilibrium situation. The correspondent number of 
phonons will show an exponential decrease: 

n (t) = ni e-fo^ dt (1) 
where rij is the number of phonons at the moment when charge equilibration begins. Fusion 



experiments seem to indicate that this will happen once the decision for fusion is taken []T3 
This means that when the system passes from the strong absorption configuration through 
the barrier the neck is large enough to allow for an isospin collective motion. Therefore we 
can assume that the configuration reached soon after, when the interdistance between the 



center of mass of the two nuclei is around the sum of their radii, represents the starting point 
of the collective pre-equilibrium GDR mode. It is important to stress that at this moment 
we have already reached a quite noticeable density overlap and therefore the considered 
configuration is not corresponding to the naive picture of the two nuclei as two touching 
rigid spheres. This will become clear from density contour plots and collective behaviours 
we will show later in microscopic dynamical simulations. 

In order to estimate rij we consider a harmonic oscillator description for the GDR: 

H Gm(t) = m + ^ (() , (2) 

where M = ^z+N) m * s ^ ne re duced mass of the neutron-proton relative motion, m = 935MeV 
being the nucleon mass, N = N1+N2, (Z = Z1 + Z2), the total number of neutrons (protons). 
Here II denotes the relative momentum: 

n = ^<§-£) (3) 

with P p (P„) the center of mass momentum for protons (neutrons), while X = R p — R n 
is the distance between the centers of mass of the two components. The initial number of 
phonons is defined as: 

H GDR (0) 1 



n, 



hw(Q) 2h 2 



M(hu(0))X 2 (0) (4) 



We use the above approximation in order to get a closed form directly related to the initial 
charge asymmetry since X(0) is given by: 

In Eq.(H) we are neglecting the kinetic part contribution. We expect this approximation not 
to work at higher incident beam energies. This contribution will be thoroughly accounted 
for from dynamical simulations, as discussed later. 



We can consider for the initial elongated shape a phonon energy pO 



MO) = IT^Egdr (6) 
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where R = r v4 1//3 , (r = 1.2), is the equilibrium radius of the CN and Egdr — 78v4 -1//3 is 
the energy of a phonon built on CN, corresponding to the centroid of the GDR spectrum 
in medium-heavy nuclei. 

In conclusion we can finally express the mean number of phonons at the time when a 
CN is formed as: 

n^DR = n(t CN ) « 7^2 MR ( R i + R ^§^^ ~ y/ EgdrGX V (-/W*Civ) 

1 A (A{ /3 + A 1 ^) Z\Z\ f N x N 2 
= 1-4 j — (- - -) exp HW^) (7) 

where p ave is an average value of the spreading width: 

1 rtcN 

Pave = / H(t)dt (8) 



tcN ^0 

We observe that v!q DR will depend critically on fusion dynamics (through the time scale for 
Compound Nucleus formation) as well as on properties of the spreading width p(t). 

After the time tcN we have to switch to the approach introduced in || for the CN decay 
phase. The number of GDR phonons is still time-dependent since charge equilibration is 
still going on, but now we can have some feeding from the CN heat bath with a width 
Ff eed = h\. The result is [§] (now the origin of time is tc7v) : 
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n GDR (t) = - 1 1 - (1 - nP DR %) exp (-/it) I (9) 

which asymptotically leads to the statistical value: 

A p(E* - Egdr) _ , „ /rr , / in s 

Jl = — ZzFj — exp ^- e gdr/t), (10) 

where p(E) is the level density, E* is the CN excitation energy and T the corresponding 
temperature. From this discussion we see that the quantity Eq.@ can be directly used as 
the dipole strength in the evaluation of 7-decay in the CN evaporation cascade. From the 
detailed balance we have a statistical 7 emission rate: 



R^^EryJ 



p{E* - E 7 ) a^ s (E 7 ) E* _p(E*-E GDR )r p{E* - £ 7 ) a abs {E^) E* 



p(E*) 3 (vrftc) 2 p{E*) Lp(E*-E GDR ) 3 (irhc) 



:n) 



where <7 a fc s (-E 7 ) is the absorption 7-dipole cross section. This leads, in presence of a pre- 
equilibrium contribution, to a time dependent quantity: 

R,(E„ t) - n GDR {t) [ p{e *_ Egdr) 3 {7rhc)2 \ (12) 

with n GDR (t) given by Eq.(|). 

The final yield of photons emitted can be then evaluated using a time-dependent evapo- 
ration cascade procedure, see ref. |ZT|. From Eq.s (|l2|) and @, once the CN initial excitation 



energy (temperature) is fixed, the results will be critically dependent on riQ , DR and Y^(T). 

To show how the method is working and the sensitivity to the above parameters in Fig.l 
we report some results for the systems studied by Flibotte et al 0]. The curves represent 
the "subtracted" photon spectra (i.e. normalized to a statistical emission without GDR): in 
(a) the data for the two systems, Ca + Mo (charge asymmetric, dashed) and S + Pd (more 
symmetric, solid), and in (b), (c), (d) the results obtained from 10 5 runs of a Monte Carlo 
evaporation cascade l^^l] with the same initial conditions for the CN and varying only 



n C!DR an d r^(T). The best results (Fig. lb) are obtained with the values Uq DR — 0.14 and 
= 8MeV. Now while it is reasonable to have a GDR damping width larger than in the 
ground state of the same compound nucleus 140 Sm {T^ = 4.8MeV used in Fig.lc) the choice 
of 

n GDR requires a deeper dynamical justification. This will be discussed later in a detailed 
connection to the fusion dynamics that will allow some interesting predictions on the beam 
energy dependence of this GDR extra strength. 

III. DYNAMICAL DIPOLE SIGNATURE ON ANGULAR DISTRIBUTIONS 

We note that the pre-equilibrium contribution has some special dynamical features, being 
associated with a dipole component on the reaction plane, actually along the symmetry axis 
of the dinuclear system. If we call this as the z— component, we expect no dipole strength 
present on the x, y components at the time of the CN formation. The corresponding number 
of phonons will also reach asymptotically the statistical value Eq.fliTiD (times 1/3 on average) 
but starting from zero at t = tcN [13 • 



We have then time dependent 7-dipole emission rates different for each component in 
the intrinsic reference system of the dinucleus: 
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- Wz{t) 3 (TTftc) 2 [L6) 



Ry^Ej, t) = Ry, y (E y , t) = — (1 - exp (-/it)) 



A 



<Jabs{E~t) E. 



\2 



3/i -I 3 (7rftc) 

" xA) 3 (vrftc) 2 l j 

with n^Q DR given by Eq.fl7|). 

We remark the physical meaning of the time dependent weights W ZjX>y (t) defined before, 
they represent the number of dipole phonons on each component at each time step. In Fig. 2 
we report a typical behaviour of the weights. We have chosen the parameters corresponding 
to the Ca + Mo system of ref. ||, see before. 

In particular we note the extreme values at t — (CN formation time) and asymptoti- 
cally t — > 00 (full charge equilibration): 

WM=n < S ) DR , W z (oc) = ±- (15) 
W XfV (0)=0 , ^(00) = ^-, (16) 

i.e. the GDR phonons are initially only on the z— component and finally uniformly dis- 
tributed on the three axes. This property will have important observable effects on the 
7— angular distributions. 



The angular distribution of the emitted photons can be written as p3J3J25.26 



W(6, E y ) = W (l + a 2 (E 1 )P 2 (6)) (17) 

where is a Legendre Polynomial in the polar angle 9 between the direction of the 

emitted 7— ray and the beam axis. With our choice of the oscillation axes we have 
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(x, z)— stretched and y— unstretched transitions |24[] and a time-dependent anisotropy pa- 
rameter: 

a ^ E ^~ w x + w y + w z ' (18) 



with the time dependent weights W ZyXyVt (t) given by Eq.s ( |T3"||H| ). The z— extra strength 
contribution should therefore lead to clear negative values around the GDR energy also in 
absence of an intrinsic deformation of the fused nucleus. We remark that all the previous 
equations can be easily extended to the use of different Lorentzians for the various axis 
contributions to a a b s , to take into account the presence of deformations. The dynamical 
dipole effect will enhance the corresponding characteristic anisotropies in the 7— angular 
correlations. 

IV. THE INTERPLAY BETWEEN FUSION DYNAMICS AND 
PRE-EQUILIBRIUM GDR 

In this section we present a fully microscopic analysis of the fusion reaction for some 
systems of interest, to look in particular at the pre-equilibrium GDR population. The calcu- 
lations are performed in the framework of the BNV transport equation which incorporates 
in a self-consistent way the mean field and two-body collisions dynamics p7| . The numerical 



accuracy has been largely improved in order to have a good description also of low energy 
fusion reactions, see discussion in |I8|,26|. Moreover in order to further reduce numerical 



fluctuations the results presented here are obtained from an average over several events. For 



the mean field we have adopted a Skyrme-like parametrization, as described in ||18|| , which 
well reproduces Nuclear Matter saturation properties (equilibrium density, binding energy, 
compressibility and symmetry energy). In the collision integral we use in medium reduced 



nucleon-nucleon cross sections, isospin as well as energy and angular dependent |28 

We have focussed our attention on the fusion reaction induced by 16 (N/Z = 1) on 98 Mo 
(N/Z = 1.33) (see J7|). We have studied in a relative large range of incident beam energy, at 
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AMeV/n, SMeV/n, lAMeV/n and 20 MeV/n, the evolution of the composite system along 
the fusion path. We will consider results that correspond to an impact parameter b = as 
well representative of the features we are looking for. 

As discussed in the previous section an essential point is to evaluate the starting time 
of phase II, when charge equilibration and the relative collective dipole are present in the 
dinuclear system. This time can be directly extracted from the collision simulation just 
comparing the evolution of the distance between Projectile- Target Centers of Mass, dpT{t), 
and of the distance between Neutron-Proton Centers of Mass, X(t), reported in Fig. 3. The 
two quantities are roughly proportional up to the time of the onset of a collective dipole 
response of the dinuclear system. At this time step, chosen as t = in Fig.3, the quantity 
X(t) clearly shows an acceleration due to the symmetry term of the dinuclear mean field. 
We remark that for all energies the value of the dpx distance corresponding to the bending 
of the "dipole" variable X(t) is very close to the sum of the two radii, in agreement with 
the discussion of the previous section. 

In all the following plots we will choose as t = such starting time of the dinuclear 
dynamics, roughly corresponding to the "touching" configuration. 

At each time step we can evaluate the dinuclear dipole moment in coordinate and mo- 
mentum space: 

DR(t) = ^X(t) (19) 
DK(t) = ^ (20) 

The results are shown in Fig. 4. The out of phase behaviour of the two dipoles, clear signature 
of a collective dinuclear response, is indeed starting just after t — (touching configuration) 
practically for all energies. We remark a smaller initial amplitude of the oscillation at the 
lowest (AAMeV) and highest (20 AMeV) energies, as well as an increase of the damping 
with the beam energy. In conclusion it clearly appears an optimum range of energies for 
the observation of the dynamical dipole effect, well above the Coulomb threshold (about 
AAMeV in the O + Mo system) but also well below the Fermi energy domain. This point 
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will be further analysed in the following. 

Fig. 5 shows phase-space trajectories of the GDR, i.e. the time evolution of the DK — DR 
correlation. The spiral curves are nicely revealing an out of phase behaviour in presence of 
some damping. The spiralling trend is starting at the touching configuration (t = 0, left 
points of the curves), maybe with some delay in the lowest energy case. The centre is 
reached when charge equilibration is achieved. From this time onwards the GDR mode will 
be only of statistical nature. The spiral shows a faster collapse to the central region at high 
beam energy (Fig.5d) since we have a larger damping of the dynamical collective motion, as 
already seen in Fig. 4. 

The amplitude reduction at 4AM eV and 20AM eV is also evident: the nature of the 
effect is however different in the two cases. At low energy, just above the Coulomb barrier, 
the neck is formed quite slowly (see Fig. 6) and so the collective dipole response of the 
dinuclear system is actually starting sometime after the touching configuration, when some 
charge equilibration has already taken place. At 20AMeV the fusion dynamics is fast but 
now we can have some prompt particle emission that will reduce the dipole moment of the 
dinucleus. Just by chance we see a quite good spiralling behaviour at 8AMeV, where the 
experiment has been performed by Cinausero et al. |7[] with a very clear evidence of the 
dynamical dipole enhancement. 

Now we can dynamically evaluate the number of pre-equilibrium GDR phonons tiq DR 
present at the time of the CN formation: as shown in the previous section this quantity 
is essentially ruling the extra emission of GDR 7-rays. At each time step the average 
number of phonons n(t) can be determined from the simulations just dividing the total 
oscillator energy, as defined by Eq.(|), by the one phonon energy at that instant, hu(t), i.e. 
corresponding to the deformation at the time t. The results are presented in Fig. 7. The 
arrows show the initial value n(0) predicted from Eq.(|j), i.e. without dynamical effects. 
The dashed curve in all four pictures is the analytical expectation of Eq.(|I]) using a constant 
spreading width value, = 5MeV. We remark a quite good overall agreement at low beam 
energies 4 and 8MeV/n, (a) and (b). With increasing energy we clearly see in the number 
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of phonons obtained from the simulations a faster damping which quickly compensates the 
higher starting point. This is particularly evident at 20MeV/n (Fig.7d) where we have 
to consider a quite larger value for the spreading width joint to a sizeable kinetic term 
contribution to the initial number of phonons n(t = 0). 

This time evolution of the number of dipole phonons in the dinuclear system will strongly 
influence the quantity we are looking for, tiq dr , given by the value of n(t) at the time of 
Compound Nucleus formation tcN, Eq.(0), the latter being fixed from equilibrium condi- 
tions in coordinate and momentum space. In Fig. 8 we plot the time evolution of the mass 
quadrupole moment in r-space at the four incident energies. Generally we remark two dis- 
tinct behaviours in the evolution: a faster decay until the first minimum is reached, followed 
by a smoother trend accompanied by small oscillations. During the first stage the conver- 
sion of relative motion into heat takes place and therefore, in the limits of some more time 
needed for a final thermal equilibration, the moment when the first minimum is attained will 
provide us the Compound Nucleus formation time. Our evaluation is in agreement with the 
information obtained by looking at the quadrupole in p-space (see Fig.9), which is a good 
probe of the energy relaxation time. We can extract tcN values of 120 fm/c at 4AMeV, 
80 fm/c at 8AMeV, and around 50 fm/c at the higher energies lAAMeV and 20AMeV. 
From Fig. 7 we get the corresponding n^Q DR = n(tcN) parameters: 0.01 at 4AMeV going 
up to 0.07 at 8AMeV and finally reducing again down to 0.05 at lAAMeV and to 0.04 at 
20AMeV. 

To test the correctness of our procedure we have performed a complete evaporation 
cascade calculation using the above extra dipole strength for the system O + Mo at 8AMeV. 
The "subtracted" spectra (see definition in section II) are shown in Fig. 10, from 2 * 10 5 
Montecarlo events, compared with experimental data from ref. [[7]. The agreement is quite 
good. 

Since the time scale of fusion dynamics is playing an important role on the amount of 
extra dipole strength present at the time of CN formation we expect to see some dependence 
of the effect also on the mass symmetry in the entrance channel, which is strongly affecting 
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the fusion time, see |21| and ref.s therein. In order to check this point we have studied the 
system 50 Cr (N/Z = 1.08) + m Ni (N/Z = 1.28) with almost the same charge asymmetry 
of O + Mo but much more symmetric in mass. We have considered central collisions at 
the beam energies 3.5 and 5.5AMeV corresponding to the same cm energies above the 
coulomb barrier of the first two cases of the O + Mo system. The relative fusion paths 
are shown in Fig.ll for 100/m/c starting from the touching configuration, to be compared 
with the equivalent Fig. 6 of the O + Mo. At both energies it is quite evident the longer 
shape equilibration time. In Fig.l2a,b (to compare with Fig.3a,b for O + Mo) we see how 
the dinuclear formation time (t = in our convention) can be again well deduced from the 
simulations. In Fig. 12c, d we show the time evolution of the correlation between momentum 
and space dipoles. If we compare with the analogous Fig.5a,b of the mass asymmetric 
O + Mo case we clearly observe a quite larger damping of the collective dipole mode in the 
dinucleus, due to the larger excitation energy available in the more symmetric system and 
the related increase of the GDR spreading width [f26],|3T]. The slower shape equilibration is 
quantitatively analysed in Fig.l3a,b (see the corresponding Fig.8a,b for O + Mo) and finally 
the time evolution of the number of dipole phonons is presented in Fig. 13c, d, to compare 
with the equivalent Fig.7a,b for the mass asymmetric case. We remark that in spite of the 
larger initial value n(t = 0), as expected from Eq.(£|) since Ri + R2 is larger in case of mass 
symmetry, we finally get a much smaller value of tiq DR = n(tcN) f° r two main reasons, the 
larger damping and the longer CN formation time. 

We have finally checked the method also for the ref. |J system, Ca + Mo at 4AM eV, 
where the best results are obtained with a tiq DR = 0.14, see Fig.l. In Fig. 14 we present 
the time evolution of the relevant quantities, mass quadrupole moment, dipole moments in 
r- and p- space and average number of phonons. With respect to the O + Mo case at the 
same energy we see a larger initial dipole strength but with a stronger damping and a slower 
fusion process. We have a tcN — 120/ra/c and a corresponding riQ DR ranging between 0.12 
and 0.15, as needed in the fitting procedure of Fig.l. 
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V. SUMMARY AND CONCLUSIONS 



We have shown that important information on the early stage of the fusion path can be 
obtained studying charge-asymmetric reactions. In fact, in this case it is possible to reveal 
a "direct" dipole oscillations (the dynamical dipole), related to the charge equilibration 
dynamics, which leads to some extra GDR strength in the statistical decay of the fused 
system. 

The dynamical nature of such pre-equilibrium contribution, i.e. corresponding to dipole 
oscillations on the reaction plane, is expected to show up in an anisotropic 7-emission, also 
in absence of deformations in the compound nucleus. 

In this work we have investigated how N/ Z asymmetry and fusion dynamics are affecting 
the properties of the pre-equilibrium GDR remnants. We have shown the existence of a 
GDR collective mode, during the formation phase of the compound nucleus, i.e. built on 
nonequilibrium states. This has allowed us to express the mean number of dipole phonons 
present in the fused nucleus in terms of the initial isospin asymmetry, CN formation time 
and GDR spreading width. We can make quite accurate predictions on the optimal choice 
of the reaction partners, in particular on the interplay between charge and mass asymmetry 
in the entrance channel. 

Particularly interesting is to follow the energy dependence of the effect. We expect to 
see a rise and fall of the dynamical dipole contribution. 

At low energies, just above the Coulomb barrier, the effect is strongly reduced for two 
combined reasons: 

• i) A delay in the dinucleus formation (slow neck dynamics) and relative collective 
response. The pre-equilibrium dipole oscillation of the composite system will have an 
initial smaller amplitude since some charge equilibration has already taken place. 

• ii) A longer Compound Nucleus formation time is decreasing the average number of 
phonons present in the fused system due to the GDR damping. 
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At higher energies, close to the Fermi energy domain, the effect is again reduced on 
average for two main reasons: 

• i) We have a relevant pre-equilibrium particle emission (i.e. incomplete fusion events) 
with a direct reduction of the initial charge asymmetry. 

• ii) In the fusion processes we have a large excitation energy deposited in the composite 
system with a related increase of the GDR spreading width leading to a fast decrease 
of the number of dipole phonons. 

Mass symmetry in the entrance channel is also strongly affecting the dynamical dipole 
extra-strength. For the same charge asymmetry and the same incident energy above the 
Coulomb barrier more mass symmetric systems are expected to show a more reduced pre- 
equilibrium GDR strength for two main reasons, a longer compound nucleus formation time 
and a larger GDR spreading width due to the larger excitation energy that can be reached 
in the dinuclear configuration. These two effects are clearly compensating the larger value 
of the dipole moment at the touching time. An increase of the beam energy will not help 
for the reasons listed above, pre-equilibrium and larger damping of the collective mode. In 
conclusion mass-symmetric entrance channels are strongly quenching the dynamical dipole 
effect at all energies. 

Moreover a new dynamical feature is appearing at beam energies above 20 AM eV, some 
large monopole oscillations in the entrance channel dynamics |29|-pT| . The dynamical dipole 



strength is reduced and more fragmented, see the simulations in ref. ||31|| , although still 
present as confirmed from very recent data pT| , [12| . The study of the dynamical dipole effect 
in the Fermi energy domain, although experimentally quite difficult, would therefore bring 
new independent information on the spreading width of hot GDR and in particular on the 
coupling to an expanding collective mode 

Finally, for charge asymmetric reactions, the enhanced dipole emission could be an in- 
teresting cooling mechanism to favour the fusion of very heavy nuclear systems. 
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From the interplay between charge and shape equilibration time-scales we can also sug- 
gest new experiments to study the dipole propagation in excited nuclei. The use of radioac- 
tive beams will enhance the possibility of such observations. 
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FIGURE CAPTIONS 



Fig. 1 Subtracted 7 spectra from a 140 Sm CN formed at E* = 71MeV in the charge asym- 
metric (Ca + Mo - dashed lines) and symmetric (S + Pd - solid lines) entrance channel: 
a) Experiment ||; b), c), d), simulations (see text) with (riQ DR , T^) respectively equal 
to (0.14, 8MeV), (0.6,8MeV), (0.14, L8MeV). 10 5 Montecarlo events. 

Fig. 2 Time evolution of the number of phonons on the different intrinsic axes for a 
140 Sm CN formed at E* = 71MeV in the charge asymmetric Ca + Mo entrance 
channel. The used parameters are Egdr = !5MeV, T = 2MeV, n^Q DR = 0.14 and 
= 8MeV. 

Fig. 3 Time evolution of the distance between Projectile- Target Centers of Mass dpr(t) 
(solid lines, left scales in fm) and of the distance between Neutron-Proton Centers of 
Mass X(t) (dashed lines, right scales in fm) for the O + Mo system at beam energies: 
a) AAMeV, b) 8AMeV, c) UAMeV, d) 20AMeV. The t = choice is discussed in 
the text. 

Fig. 4 Time evolution of the dipole moment in r-space, DR (solid lines, left scales in fm) 
and in j>-space, DK (dashed lines, right scales in fm' 1 ) for the O + Mo reaction. 
t = corresponds to the touching configuration (see text). The labels correspond to 
the same energies of Fig. 3. 

Fig. 5 Phase space trajectory of the entrance channel dipole for the O + Mo reaction. The 
labels correspond to the same energies of Fig. 3. 

Fig. 6 Density plots of the neck dynamics for the O + Mo system at the two energies 
4AMeV and 8AMeV. 

Fig. 7 Time evolution of the number of dipole phonons for the O + Mo reaction (solid 
lines). The labels correspond to the same energies of Fig. 3. The dashed line is the 
reference curve discussed in the text, with a constant spreading width. 
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Fig. 8 Time evolution of the mass quadrupole moment for the O + Mo reaction (in a.u.). 
The labels correspond to the same energies of Fig.3. 

Fig. 9 Time evolution of the quadrupole moment in momentum space for the O + Mo 
reaction (in a.u.). The labels correspond to the same energies of Fig.3. 

Fig. 10 Subtracted 7 spectra from a U4 Sn CN formed at E* = 108MeV in the charge 
asymmetric (O + Mo - dashed lines) and symmetric (Ti + Ni - solid lines) entrance 
channel: a) Experiment @; b) simulations (see text) with (tiq DR , F^) respectively 
equal to (0.07, 7MeV). 2 * 10 5 Montecarlo events. 

Fig. 11 Density plots of the neck dynamics for the Cr + Ni system at the two energies 
3.5AMeV and 5.5AMeV. 

Fig. 12 Collision of the more mass symmetric Cr + Ni system at the two beam energies 
3.5 and 5.5AMeV . (a) and (b): quantities like in Fig.3. (c) and (d): quantities like in 
Fig.5. 

Fig. 13 Collision of the more mass symmetric Cr + Ni system at the two beam energies 
3.5 and 5.5AMeV. (a) and (b): quantities like in Fig. 8. (c) and (d): quantities like in 
Fig.7. 

Fig. 14 Time evolution of mass quadrupole moment (a), r-dipole moment (b), j>-dipole 
moment (c) and average number of dipole phonons (d) for the reaction Ca + Mo at 
AAMeV. 
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Yields (a.u.) 



o m. to CO ^. Oi 




O >->■ tO CO ^ Oi 



